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ABSTRACT 

We use the Millennium Simulation series to investigate the relation between the ac- 
cretion history and mass profile of cold dark matter halos. We find that the mean 
inner density within the scale radius, r_2 (where the halo density profile has isother- 
mal slope), is directly proportional to the critical density of the Universe at the time 
when the virial mass of the main progenitor equals the mass enclosed within r_2- 
Scaled to these characteristic values of mass and density, the average mass accretion 
history, expressed in terms of the critical density of the Universe, M (p C rit{z)) , resem- 
bles that of the enclosed density profile, M((p)), at z = 0. Both follow closely the 
NFW profile, which suggests that the similarity of halo mass profiles originates from 
the mass-independence of halo accretion histories. Support for this interpretation is 
provided by outlier halos whose accretion histories deviate from the NFW shape; their 
mass profiles show correlated deviations from NFW and are better approximated by 
Einasto profiles. Fitting both M((p)) and M(p cr it) with either NFW or Einasto pro- 
files yield concentration and shape parameters that are correlated, confirming and 
extending earlier work that has linked the concentration of a halo with its accretion 
history. These correlations also confirm that halo structure is insensitive to initial con- 
ditions: only halos whose accretion histories differ greatly from the NFW shape show 
noticeable deviations from NFW in their mass profiles. As a result, the NFW profile 
provides acceptable fits to hot dark matter halos, which do not form hierarchically, 
and for fluctuation power spectra other than CDM. Our findings, however, predict a 
subtle but systematic dependence of mass profile shape on accretion history which, if 
confirmed, would provide strong support for the link between accretion history and 
halo structure we propose here. 
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1 INTRODUCTION 

Numerical simulations have shown that the equilibrium 
structure of cold dark matter (CDM) halos is approximately 
self-similar. Spherically averaged density profiles, in partic- 
ular, are well approximated by scaling a simple formula pro- 
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posed by iNavarro et al.l (| 19951 . Il99rj . hereafter NFW). The 
NFW profile has fixed shape, and is characterized by a log- 
arithmic slope that steepens gradually from the center out- 
wards. As such, it may be fully specified by just two pa- 
rameters, usually chosen to be either the virial radius and a 
characteristic density or, equivalently, the halo virial mass 
and a concentration parameter. (See Sec. 13. f I for a summary 
of refevant formulae and definitions.) 
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The gently-varying slope 

of the NFW profile confounded early theoretical expecta- 
tions, whj£h_Jiad_ejiyjsioMd^^ 

fe.g.,lFillmore fc Goldreich]|l984l; iHoffman fc Shahamll 19851 ; 



h alos selected from th e three Millennium Simulations, MS - 
I (|Springel et al.ll2005l), MS-II dBoylan-Kolchin et alj|2009h . 



lQuinn~t al.l 19861 ; Crone et al.lll994 ). and has motivated 
a number of proposals to explain its origin (see, for a 
recent review, iFrenk fc White! |2012| ). Most attempts have 
ta ken as guidance the se condary-infall model first proposed 
by iGunn fc GottJ (|l972T ) , complemente d by various conjec- 



tures about the role of merger s (e.g., Salvador-Sole et all 



Il998l ). dynamical fricti on (e.g., Nusser fc Shetbl 19991 ). an 



gular momentu m (e.g.. I Williams et all 120041'). "or adiabati c 
invariants (e.g.. lAvila- Reese et al.l 19981 ; Dalai et all |2010T ) . 
or else have argued that entropy generation during virial- 
izati on might be behind the halo structural similarity (s ee, 
e.g.. lTavlor fc Navarrdl200ll ; | Pontzen fc Governatoll2013l ). 

No general consensus has yet emerged, however, reflect- 
ing the difficulty that all of these models face when trying to 
explain why the same NFW profile seems to fit the structure 
of halos formed t hrough hierarchical clustering regardless of 
power spectrum (|Navarro et al.|[l997l ). as well as that of hot 
dark matter halos or systems formed through mon olithic 
collapse (e.g., lHuss et~al]|l999l ; IWang fc Whitell2009l ). 

In addition, none of these models provides a thor- 
ough explanation for the redshift-dependent correlations be- 
tween mass and concentration seen in simulations, their 
scatter, or their dependence on cosmological parameters 
and power spectra. Halo concentration, which depends 
only weakl y on mass, was or iginally linked to halo col- 
lapse time (jNavarro et al.|[l997t l . but attempts to reproduce 
the simulation results with simple prescript ions based on 
that proposal have met with limited success jBullock et al.l 
20011; lEke et al.ll200ll; fNeto et al.ll2007l; iMaccio et al.ll2008l; 



Gao et al.l l2008) 



Better results have been obtained with empirical mod- 
els that relate concentration to halo mass accretion history 
and, in particular, to the time when the main halo progeni- 
tor switc hes from a period of " f ast growth" to one of "slow 
grow th" (|Wechsler et all |2002| ; IZhao et al.l 120031 ; iLu et al.1 
l2006h. The succ e ss of these models is not, however, unquali- 
fied. Zhao et ail (|2009f ) , for example, argue that halo concen- 
tration is determined at the time when the main progenitor 
first reaches 4% of the final mass, but there seems to be no 
natural justification for why concentration should be related 
to this particular, rather arbitrary time of a halo's assembly 
history. 

Further complicating matters, there is now convinc- 
ing evidence that a number of halos have density profiles 
that deviate slightly, but significantly, from the NFW pro- 
file (|Navarro et al.l 1 2004 '). Accounting for these deviations 
requires the introduction of an additional shape parameter, 
thus breaking the structural similarity of CDM halos. One 
parameterization that results in excellent fits is the Einasto 
profile, where the logarithmic slope is a simple power law 
of radius, dlnp/dlnr cx (r/r_2) a : the shape parameter, 
a, is the exponent of the power l aw. This finding has now 
been verified by additional work dMerritt et alj|2005l. I2OO6I; 



Navarro et alj|2010; Gao ct al. 2008; Havashi & White 2008 



Stadel et ai]|2009l ; iLudlow et al]l201ll ') but there is no clear 



understanding of what breaks the similarity or what deter- 
mines the value of a for a particular halo. 

We explore these issues here using a large ensemble of 



and MS-XXL (|Angulo et al.l 12012? ) . collectively referred to 
hereafter as MS. These are amongst the largest cosmologi- 
cal A r -body simulations available, and provide us with thou- 
sands of well-resolved halos spanning more than four decades 
in mass. Merger trees are available for all these simulations, 
making them an ideal dataset to explore the relation be- 
tween accretion history and mass profiles. In addition, the 
numerical homogeneity and sheer size of the volumes sur- 
veyed by the MS allow us to combine large numbers of halos 
with similar properties to smooth out statistical fluctuations 
and idiosyncrasies of individual systems that might obscure 
the general trends. Our analysis reveals a subtle but well- 
defined relation between mass profile and accretion history 
that offers valuable new clues to the origin of the structure 
of CDM halos. 

The plan of this paper is as follows. We describe the 
simulations in Sec. [2] and the analysis procedure in Sec. [3] 
We present our main results in Sec. [4] and summarize our 
main conclusions in Sec. [S] 



2 NUMERICAL SIMULATIONS 

Our analysis focuses on dark matter halos identified in the 
three Millennium Simulations. We provide here a brief sum- 
mary of these simulations and of their associated halo cata- 
logs. We refer the reader to the original papers for extensive 
details on each of the MS runs. 



2.1 The Millennium Simulation suite 

All MS runs adopt a flat, WMAP 1-normalized LCDM cos- 
mology with the following cosmological parameters: £7 m = 
0.25, Ov = 1 - fi m = 0.75, h = 0.73, n = 1 and cr 8 = 0.9. 
Here S7; is the present-day contribution of component i to 
the total matter energy density in units of the critical den- 
sity for closure, p C rit; is the rms mass fluctuation in 
8 ft -1 Mpc spheres, linearly extrapolated to z = 0; n is 
the spectral index of primordial density fluctuations, and h 
is the Hubble parameter. In addition to using the same cos- 
mological parameters, the MS runs also adopted the same 
sequence of outputs in order to facilitate comparisons be- 
tween the runs. 

MS-II follows the dark matter distribution using 2160 3 
particles of mass m p = 6.89 x 10 6 h~ x M @ in a 100/i _1 Mpc 
periodic box. MS-I has the same total particle number, but 
follows the evolution of structure in a comoving box of 
500 h" 1 Mpc on a side; each particle in MS-I is thus 125x 
more massive than in MS-II, or m p = 8.61 x 10 s ft -1 Mq. 
MS-XXL is the largest of the three simulations in both box 
size and particle number; it follows 6720 3 particles of mass 
m p = 6.17 x 10 9 /i _1 M® in a 3 /i^Gpc box. 



2.2 Halo Catalogs 

A friends-of- friends (FOF) group finder (|Davis et al.lll985h 
was run on the fly for each simulation output using a linking 
length of b — 0.2 times the mean inter-particle separation 
and a minimum particle number iV m i n = 20. The subhalo 
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Figure 1. Halo density profiles and accretion histories. Left panel: Median density profiles of MS-1I relaxed halos in the mass range 
f.24 < log M2oo/(lO lo /i -1 M0) < f.54 (corresponding to particle numbers in the range 2.5 X fO 4 < ./V200 < 5 X 10 4 ), selected according 
to their concentration (see boxes in the top panel of Fig. [2j. Densities arc shown scaled to po, the critical density at z = 0, and weighted 
by r 2 in order to enhance the dynamic range of the plot. Radii are scaled to the virial radius, r200- The best-fit Einasto profiles arc shown 
by the thin solid curves, with parameters listed in the legend. Dot-dashed curves indicate NFW profiles (whose shape is fixed in these 
units) matched at the scale radius, r_2, where the r 2 p profiles peak. Arrows indicate the half-mass radius, rj/2- Right panel: Median 
mass accretion histories (MAH) of the same set of halos chosen for the left panel. Halo accretion history is defined as the evolution of 
the mass of the main progenitor, expressed in units of the mass of the halo at z = 0. The heavy circles indicate the redshift, z_ 2, when 
the progenitor's mass equals the mass, M_2, enclosed within the scale radius at z = 0. Starred symbols indicate the half-mass formation 
redshift. 



finder SUBFIND, (|Springel et al.|[200lh was then run to iden- 
tify self-bound substructure within each FOF halo. 

Subfind dissects each FOF halo into one dominant 
structure (the main halo) and a number of subhalos that 
trace the self-bound remnants of accreted systems. We will 
focus our analysis on main halos identified at z = that 
contain at least -/V200 = 5000 particles within their virial 
radiufQ. 

Since dark matter halos are dynamical systems, tran- 
sients induced by mergers or ongoing accretion can lead to 
rapid fluctuations in the structure of a halo that are poorly 
captured with simple fitting formulae. We therefore impose 
three criteria to flag systems that are clearly out of equilib- 
rium. We consider a halo to be dynamically "relaxed" if it 
simultaneously sa tisfies all three of the following conditions 
(|Neto et al.ll2007l ): (i) / sub < 0.1, (ii) d oS < 0.07 and (hi) 
2T/\U\ < 1.35. Here / su b is the fraction of the halo's virial 
mass contributed by subhalos, d g = |r p — rcivi|/r2oo is the 
distance between the halo barycenter and the location of its 
potential minimum, expressed in units of r2oo; and 2T/\U\ is 
the virial ratio of kinetic to potential energies, measured in 
the halo rest frame. None of our conclusions are heavily af- 
fected by these restrictions. Unrelaxed systems make up only 
20% of all halos with virial mass of order 10 12 M and 25% 



We define the virial radius, r2oo, of a halo as the radius of a 
sphere centered at the potential minimum that encloses a mean 
density of 200 X p cr it- We identify all virial quantities (i.e., mea- 
sured within r20o) with a "200" subscript. 



of ~ 10 Mq halos. Only at very large halo masses, such as 
cluster-sized ~ 10 14 Mq systems, the unrelaxed frac tion ex- 
ceeds 50%. We refer the reader to lNeto et~al"1 (l2007h for fur - 
ther discussion of these criteria, and to Ludlow et alj (|2012l ) 
for a discussion of how the inclusion of out-of-equilibrium 
systems may impact the mass-concentration relation at large 
halo masses. 



3 ANALYSIS 

3.1 Fitting Formulae 

We consider two different formulae to fit halo density pro- 
files. The NFW profile is given by 



(1) 



P(r) = 8c 

pcrit (r/r s )(l + r/r s ) 2 ' 

where r B is a scale radius, p cr j t = 3H 2 /8ttG is the critical 
density, and S c is the halo dimensionless characteristic den- 
sity. These two parameters can also be expressed in terms of 
the halo virial mass, M200, and a concentration parameter, 
c = T2oo / r s , which is related to S c by 



200 



3 [ln(l + c) - c/(l + c)] 



(2) 



Note that for given mass the NFW profile has a single free 
parameter, the concentration. This profile can also be ex- 
pressed in terms of the enclosed mean density, M((p}), where 
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<p)(r) = 



M(< r) _ 200 Y(cx) 
(47r/3) r 3 ~~ a; 3 F(c) 



Pcrit , 



(3) 
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Figure 2. Mass dependence of the best- fit Einasto parameters 
for all halos in our sample at z = 0. Only relaxed halos with 
more than 5000 particles within the virial radius are considered. 
The top and bottom panels show, respectively, the concentration, 
c = r20oA-2, and shape parameter, a, as a function of halo 
virial mass. Individual points are colored according to the third 
parameter (see color bar on the right of each panel). Connected 
symbols trace the median values for each Millennium Simulation 
(see legend in the top panel); thin solid lines delineate the 25 
to 75 percentile rang e. The dashed cu rves indicate the fitting 
formulae proposed bv lGao et al. | Hooi). For clarity only 10,000 
halos per simulation are shown in this figure. Halos shown in 
grey are systems where the best-fit scale radius is smaller than 
the convergence radius; these fits are deemed unreliable and the 
corresponding halos are not included in the analysis. Grey vertical 
bars highlight three different mass bins used to explore parameter 
variations at fixed halo mass (see Sec. 14.31 and 14,41 . Small boxes 
indicate halos in each of those bins with average, highcr-than- 
average, and lower-than-average values of a (bottom panel) or of 
the concentration (top panel). 



where x = r/r2oo and Y( u) = ln(l + u) — u/(l 4- u). 

The Einasto profile l|Ein asto 1965) has an extra free 
parameter, the shape parameter a, and may be written as 



In 



(£) 




2 








a 





(4) 



The parameter r_2 marks the radius where the logarithmic 
slope of the density profile is equal to —2. The same property 
holds for the NFW scale radius, r s , and therefore, for short, 
we shall hereafter refer to the scale radius of either profile as 
r_2- Quantities measured at (or within) r_2 will be denoted 
by a "—2" subscript; e.g., {p-2) = (p)(r_2). Of course, like 
NFW, the Einasto profile may also be expressed in terms of 
its enclosed mean density profile, M((pe))- 

We note that, for given concentration, an Einasto profile 
with a « 0.18 resembles closely an NFW profile over roughly 
two decades in radius or enclosed mass. Profiles with other 
valu es of a deviate systematica lly from the NFW shape (see, 
e.g., lNavarro et al .11 20041 , |201Ch . 



3.2 Profile Fitting 

Our analysis deals primarily with the spherically averaged 
density profiles of relaxed CDM halos identified at z — 
in each MS. We construct radial profiles using 32 concen- 
tric bins, equally spaced in logr, spanning the radial range 
-2.5 < logr/raoo < 0. 

The Einasto profile has three free parameters: p-2, T-2, 
and a. These are simultaneously adjusted in order to min- 
imize its rms deviation from the binned density profiles. In 
practice, we define a figure of merit, 



IN Qui 

^ 2 = ivb^ [lnpi 



In p B (p_ 2 ;r_2;a)] 2 , 



(5) 



which is minimized to obtain the best-fitting set of param- 
eters for any given halo. In practice, the parameters p_2 
and r_2 can be expressed in a variety of equivalent forms, 
such as virial mass and concentration (M2oo,c), or the mag- 
nitude and location of the peak in the circular velocity curve 
(V max , r max ). In order to ease comparisons with previous 
work, we characterize the dark matter halo mass profile in 
terms of its virial mass Mo = M2oo(z = 0), its concentra- 
tion c = r2oo/f-2, and its Einasto "shape" parameter, a. 
The Einasto profile provides an excellent description of the 
density profile of relaxed MS halos: the median value of ip 
is just 0.073i{5"oii> where the range represents the 25th and 
75th percentiles. 

An analogous procedure is used when NFW fits need to 
be performed; in this case, the two parameters estimated by 
the fit can also be expressed as the virial mass and concen- 
tration. 

The fits are carried out over a radial range r m i n < r < 
f200- The fitting procedure yields robust estimates for p-2, 
r_2 and a, provided r m i n is chosen to be the minimum of 
either r conv or .05 x r2op. Here, r CO nv is the convergence 
radius defined bv lPower et al.l (|2003l ). where circular velocity 
profiles converge to better than ~10%. 
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Figure 3. Relation between mass profiles at z = and accretion histories for relaxed, well-resolved halos (-/V200 > 2.5 X 10 4 ) in our 
sample. Individual halos are colored by mass, according to the color bar at the top of the plot. Left panels: Mean enclosed densities 
within the radii, ri/2, and r 3 /4, containing, respectively, 25%, 50% and 75% of the virial mass, shown as a function of the (critical) 

density of the Universe at the time when the progenitor's virial mass equals the mass enclosed within each of those radii at z = 0. These 
densities are correlated, as expected if denser halos collapse earlier. However, the dependence varies with radius and is generally quite 
weak. This explains, for example, why measures of halo density (such as the concentration) correlate only poorly with the half-mass 
formation time. Medians, quartiles, and 10/90 percentiles are indicated by the box-and-whisker symbols. Right panels: As the left panels, 
but for radii equal to half, one, and two times the scale radius, r_2. The dotted line indicates direct proportionality, scaled vertically to 
best fit the data of each panel (fit parameters given in the legends). The excellent agreement between this simple scaling and the data 
implies that, expressed in units of the scale radius, the shape of the mass profile of a halo is intimately related to that of the accretion 
history of its main progenitor. 



3.3 Mass Profiles and Accretion Histories 

The left panel of Fig. [T] illustrates the role of c and a in 
describing the density profile. This figure shows the den- 
sity profile of MS-II halos selected in a narrow range of 
mass, 1.24 < log Af 2 oo/lO lo h _1 M0 < 1.54. (Densities are 
weighted by r 2 in order to enhance the dynamic range of the 
plot.) Each profile corresponds to different systems, grouped 



by concentration: green squares track the mediaijf] profile of 
halos with average concentration for that mass; blue circles 
and red triangles correspond to halos with concentration 
~ 50% higher and lower than the average, respectively (see 
boxes in the top panel of Fig. [2]). 

2 Median profiles are computed at each radius after scaling all 
individual profiles as in Fig. [T] 



6 Ludlow et al. 




Figure 4. Average mass profiles at z = and accretion histories for halos in three different mass bins (see shaded regions in the bottom 
panel of Fig.[2jl- Top left: Average mass profiles of all halos in each bin, plotted as enclosed mass (in units of M2oo)j versus inner density 
(in units of 200x the critical density). Dashed lines are best-fit NFW profiles, which have a single adjustable parameter, the concentration, 
c = r20oA— 2- Heavy symbols indicate the enclosed mass, M_2, and density, (p_2), at the scale radius of each profile. Residuals from 
the best fits are shown in the bottom inset. Top right: Same as top-left panel, but scaled to the enclosed mass, M_2, and ovcrdensity, 
(p-2), at the scale radius. Scaled in this manner, halo mass profiles all look alike and are very well approximated by an NFW profile 
(dashed curve). Bottom left: Average accretion histories of the same halos shown in the top panels. The plots show the growth of the 
virial mass of the main progenitor, normalized to the final mass at z = 0, as a function of time, expressed in terms of the critical density 
of the Universe at each redshift. The dashed curves are not fits to the data. Rather, they indicate accretion histories parameterized, as 
in the top panel, by an NFW profile in this M-p plane. The single adjustable parameter to these profiles is fully specified by the filled 
heavy symbols, which indicate M_2, chosen to match that of the mass profiles (top-left panel) and by p C rit (•£— 2)1 computed as 776 (p— 2) 
following the correlation shown in the middle panel of Fig. [3] The open heavy symbols indicate the scale mass and density of the fitted 
NFW profile. Bottom right: Same accretion histories as in the bottom-left panel, but scaled to M_2 and (p_ 2) (open heavy symbols in 
the bottom-left panel). Note the remarkable similarity in the shape of the halo mass profiles at z = and that of the accretion histories 
of their main progenitors. 
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In the scaled units of Fig. [T]the scale radius, r_2, sig- 
nals the location of the maximum of each curve, and dif- 
ferent concentrations show as shifts in the position of the 
maxima, which are indicated by large filled circles. In ad- 
dition to their different concentrations, the profiles differ as 
well in a, which increases with decreasing concentration (see 
legends in Fig. [!}. Arrows indicate the half-mass radius of 
each profile. Dot-dashed curves show NFW profiles (whose 
shape is fixed in this plot) with the same concentration as 
the best Einasto fit (solid lines). The density profile curves 
more gently than NFW for a £ 0.18 and less gradually than 
NFW for a ^ 0.18, respectively. 

The (median) mass accretion histories corresponding to 
the same sets of halos are shown in the right-hand panel of 
Fig. [1] We define the mass accretion history (MAH) of a 
halo as the evolution of the virial mass of the main pro- 
genitor, usually expressed as a function of the scale factor 
a = 1/(1 + 2:), and normalized to the present-day value, 
Mo = M2oo(z = 0). As expected, more concentrated ha- 
los accrete a larger fraction of their final mass earlier on. 
Filled stars indicate the "half-mass formation redshift" , Z1/2, 
whereas filled circles indicate z-2, the redshift when the 
mass of the main progenitor first reaches M_2, the mass 
enclosed within r_2 at z = 0. 



1.4 



4 RESULTS 

4.1 The mass-concentration-shape relations 

The top panel of Figure [2] shows the mass-concentration re- 
lation for our sample of relaxed halos at z — 0. Concentra- 
tions are estimated from Einasto fits, and are color coded by 
the shape parameter, a, as indicated by the color bar. Open 
symbols track the median concentrations as a function of 
mass. Thin solid lines trace the 25th and 75th percentiles 
of the scatter at fixed mass. Different symbols are used for 
the different MS runs, as specified in the legend. Note the 
excellent agreement in the overlapping mass range of each 
simulation, which indicates that our fitting procedure is ro- 
bust to the effects of numerical resolution. 

The bottom panel of Fig. [2] shows the mass-a relation, 
colored this time by concentration. The trend is again con- 
sistent with earlier work; the median values of a are fairly in- 
sensitive to halo mass, except at the highest masses, where it 
increases slightly. The mass-concentration-shape trends are 
consistent with earlier work; for example, the d ashed lines 
corres pond to the fitting formulae proposed by iGao et all 
|2008l ) and reproduce the overall trends very well. 

Fig. [2] illustrates an interesting point already hinted at 
in Figure [1] the shape parameter seems to correlate with 
concentration at given mass. Interestingly, halos of average 
concentration have approximately the same shape parameter 
(a ~ 0.18, i.e., quite similar to NFW), regardless of mass. 
Halos with higher-than-average concentration have smaller 
values of a and vice versa. This suggests that the same mech- 
anism responsible, at given mass, for deviations in concen- 
tration from the mean might also be behind the different 
mass profile shapes at z = parameterized by a. We ex- 
plore this possibility next. 
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Figure 5. Relation between concentration parameters obtained 
from NFW fits to the average accretion histories and mass pro- 
files shown in Fig. fJ] The dashed curve indicates the expected 
concentration-concentration dependence given the correlations 
shown in the middle-right panel of Fig. [3] the dotted line shows 
the simple relation c [M(p)] = 2 (1 + c [MAH]), which approxi- 
mates the numerical solution rather well. Note that the relation 
is rather shallow, indicating that even halos whose accretion his- 
tories differ greatly may have similar concentrations, a result con- 
sistent with the weak mass-concentration dependence reported in 
earlier work. 



4.2 Characteristic densities and assembly times 

As pointed out by iNavarro et all (119971 ) and confirmed by 
subsequent work (see, e.g., Jing 200oh . the scatter in concen- 
tration is closely related to the accretion history of a halo: 
the earlier (later) a halo is assembled the higher (lower) its 
concentration. 

This is clear from the assembly histories shown in Fig. [1] 
which illustrate as well that defining "formation time" in a 
way that correlates strongly and unequivocally with concen- 
tration is not straightforward. For example, the often-used 
half-mass formation redshift, Z1/2, varies only we akly with c, 
maki ng it an unreliable proxy for concentration (|Neto et al.l 
1 2007f ) . An ideal definition of formation time would result in 
a natural correspondence between the characteristic density 
of a halo at z = and the density of the Universe at the 
time of its assembly. 

We explore two possibilities in Fig [3] Here we show the 
mean density enclosed within various characteristic radii at 
z — versus the critical density of the Universe at the time 
when the main progenitor mass equals the mass enclosed 
within the same radii. 

The left panels correspond to radii enclosing 1/4, 1/2, 
and 3/4 of the virial mass of the halo. Dots indicate individ- 
ual halos colored by halo mass, as shown in the color bar at 
the top. Boxes and whiskers trace the 10th, 25th, 75th, and 
90th percentiles in bins of p cr it. Note the tight but rather 
weak (and non-linear) correlation between densities at these 
radii. This confirms our earlier statement that "half-mass" 
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Figure 6. Mass accretion histories, M(p CI ^ t ), corresponding to 
Einasto profiles, compared with NFW. Note that NFW resembles 
closely an Einasto profile with a ~ 0.18 or so. Larger or smaller 
values of a correspond to halos that have been assembled more 
or less rapidly than the NFW curve, respectively. Residuals from 
NFW are shown in the bottom panel. 

formation times are unreliable indicators of halo character- 
istic density: halos with very different Z\i% may nevertheless 
have similar concentrations. 

The right-hand panels of Fig. [3] show the same den- 
sity correlations, but measured at various multiples of r_2, 
the scale radius of the mass profile at z = 0. The mid- 
dle panel shows that the mean density within r_2, (p-2) = 
-M-2/(47r/3)ri2, is directly proportional to the critical den- 
sity of the Universe at the time when the virial mass of the 
main progenitor equals Af_2. Intriguingly, this is also true at 
r-2/2 (top right panel) and at 2 x r-2 (bottom right panel), 
although with different proportionality constants (listed in 
the figure legends). 

This means that there is an intimate relation between 
the mass profile of a halo and the shape of its mass ac- 
cretion history, in the sense that, once the scale radius is 
specified, the MAH can be reconstructed from the mass 
profile, and vice versa. Since mass profiles are nearly self- 
similar when scaled to r_2, this implies that accretion his- 
tories must also be approximately self-similar when scaled 
appropriatel y. The MAH sel f- simila rity has been previously 
discussed by Ivan den Boschl (|2002h , but its relation to the 
shape of the mass profile, as highlighted here, has so far not 
been recognized. 

4.3 NFW accretion histories and mass profiles 

We explore further the relation between MAH and mass pro- 
file by casting both in a way that simplifies their comparison, 
i.e., in terms of mass versus density. In the case of the mass 
profile, this is just the enclosed mass-mean inner density re- 
lation, M{{p)) (see Sec. 13 . 1 P> . For the MAH, this reduces to 
expressing the virial mass of the main progenitor in terms 



of the critical density, rather than the redshift, M(p cr it(^)). 
In what follows, we shall scale all masses to the virial mass 
of the halo at z = 0, Mo; p C rit(z) to the value at present, po; 
and (p) to 200 po. 

The top-left panel of Fig. [Jshows, in these scaled units, 
the average M({p)) profile for halos in three different nar- 
row mass bins (indicated by grey vertical bars in the bot- 
tom panel of Fig. [2]). These mean profiles are computed by 
averaging halo masses, for given (p), after scaling all indi- 
vidual halos as indicated above. As expected, each profile is 
well fit by an NFW profile where the concentration increases 
gradually with decreasing mass. The heavy symbols on each 
profile indicate the value of M-2 and {p-i}. The top-right 
panel shows the same data, but scaled to these character- 
istic masses and densities. Clearly the three profiles follow 
closely the same NFW shape, which is fixed in these units. 

The corresponding MAHs, computed as above by av- 
eraging accretion histories of scaled individual halos, are 
shown in the bottom-left panel of Fig. [4] The heavy symbols 
on each profile again indicate the value of M-2 (as in the 
above panel), as well as p C rit(z-2) = 776 (p-2), computed 
using the relation shown in the middle-right panel of Fig. [3] 

In these scaled units, a single point can be used to spec- 
ify the "concentration" of an NFW profile, which is shown by 
the dashed curves. Interestingly, these provide excellent de- 
scriptions of the MAHs: rescaled to their own characteristic 
density and mass they all look alike and also follow closely 
the NFW shape (bottom-right panel of Fig. Q. The mass 
accretion histories and mass profiles of CDM halos are not 
only nearly self-similar: they both have similar shapes that 
may be approximated very well by the NFW profile. 

This implies that the concentration of the mass pro- 
file just reflects the "concentration" of the MAH. Indeed, 
assuming that the NFW shape holds for both, the relation 
p C iit{z-2) = 776 {p-2) delineates a unique relation between 
the two concentrations, which is shown as a dashed line in 
Fig. [5] The three symbols in the same figure correspond to 
the three average profiles and MAHs shown in Fig. 0] and 
clearly follow the same relation. 

Note that this is consistent with earlier claims that 
halo concentration is linked to the time when halo 
growth switches from a fast- to a s low-accretion phase 
(Wechsler et al.ll2002i ; Tzhao et al.ll2003l l. In our interpreta- 
tion, since both the MAH and the mass profile follow the 
same NFW shape the scale radius of one tracks that of the 
other: the "curvature" of the MAH is therefore reflected in 
that of the mass profile. Note as well that the relation shown 
in Fig. [S] is rather weak; in other words, even large changes 
in the MAH map onto a small range of concentrations in the 
mass profiles. This is at the root of the weak correlation be- 
tween co ncentration and v irial mass reported in earlier work 
(see, e.g-. lNeto et~ai]|2007| ). 



4.4 Einasto accretion histories and mass profiles 

The striking similarity between the shapes of the MAH and 
mass profile discussed above suggests an explanation for why 
halos that are outliers in the mass-concentration relation 
tend to have mass profiles that differ more significantly from 
NFW (i.e., they have a parameters that differ from 0.18, see 
Fig- El) • In this interpretation, outliers in M200-C have MAH 



The Mass Profile and Accretion History of Cold Dark Matter Halos 9 



log M 2M =[1.24, 1.54] : 

a = 0.248, c = 8.0 

a = 0.171, c=10.7 ^ 

a = 0.094, c=15.9 \ 

" Best-fit NFW \ 

\ 


. log M 200 =[3.33, 3.63] _ 

- a = 0.253, c = 5.5 
! ex=0,175, c = 6.4 

- a = 0.097, c = 7.9 % 

v . . . . 


log M 2M =[4.19, 4.49] I 

- a = 0.275, c = 4.5^^ 

! a = 0.196, c = 5.0 % 

- a=0.H9, c = 5.1 


■>:!!!::::■::;:!: 

: Residual from best-fit NFW 







0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

log <p>/(200x Po ) log <p>/(200x Po ) log <p>/(200x Po ) 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

log p crlt (z)/p log p crit (z)/p log p crit (z)/p 



Figure 7. As Fig. [4] but for halos with highcr-than-average (blue), average (green), or lowcr-than-average (red) values of the Einasto 
parameter a (see boxes in the bottom panel of Fig. [2}. Left, middle, and right panels correspond to each of the three mass bins, as 
indicated in the legends. Top panels: Average mass profiles compared with the best-fit NFW profile for all halos of the same mass (see top 
left panel of Fig. [4j . Residuals from that profile are shown at the bottom of each panel. Note the similarity between the residual curves 
of similar color at all masses. Different values of a imply different profile shapes, and deviate systematically from NFW. Bottom panels: 
Average mass accretion histories corresponding to the same halos as in the top panels. The dashed curves indicate the average "NFW 
accretion histories" for each mass bin, as shown in the bottom-left panel of Fig. [4] Residuals from this average history are shown in the 
bottom inset of each panel. Note the similarity between the shape of the residual curves of similar colors in all panels. This indicates 
that the mass accretion history is intimately linked to the mass profile at z = 0. Halos that, at z = 0, have mass profiles that deviate 
from NFW in a particular way have accretion histories that deviate from the NFW shape in a similar way. 



shapes that differ systematically from the mean, NFW-like 
shape. 

In order to test this, we may use the Einasto formula to 
fit both MAH and mass profiles. Fig. [6] shows Einasto M- 
p profiles for various values of ot, and compares them with 
an NFW profile of the same concentration. (This figure uses 
the same scalings as Fig. [4]) As stated earlier, over the range 
of mass and density plotted here the NFW profile is essen- 



tially indistinguishable from an a = 0.18 Einasto profile, 
but systematic deviations become apparent for other values 
of a. Interpreting Fig. [6] as a mass accretion history, we see 
that a > 0.18 corresponds to halos that are assembled more 
rapidly than expected from the NFW shape. The opposite 
holds for a < 0.18. This behaviour is clearly seen in the 
residuals from NFW, which are shown in the bottom inset 
of the figure. 
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Figure 8. Concentrations and shape parameters of Einasto profiles fitted to either accretion histories or mass profiles at z = 0. Heavy 
symbols correspond to well-resolved halos grouped according to the c and a parameters of their mass profile (see details in the text). 
Grey dots correspond to individual halos in the same three mass bins chosen fin Fig. [7J The left panel shows that the shape of the mass 
accretion history and that of the mass profile are correlated. The panel on the right is analogous to Fig.[5]and shows that the same applies 
to the concentrations. In this case, the relation depends on the value of a, as shown by the colored lines labelled in the legend. The heavy 
symbols are of the same type as in the left panel, but colored by a (see inset). Note that the correlations are relatively shallow, implying 
that even large departures from NFW-like mass accretion histories lead only to minor deviations from NFW in the mass profiles. 



The top panels of Fig. [7] show the average M((p)) pro- 
files of halos in three different narrow mass bins chosen to 
have different values of a. These are halos whose Einasto 
parameters fall in the boxes drawn in the bottom panel of 
Fig. [2] The best-fit NFW profile for each mass bin (as in 
Fig. [4| is indicated by a dashed curve in each panel. Devia- 
tions from the NFW curve are shown in the residuals panels. 
As expected, the residuals have different shapes depending 
on the value of their shape parameter a. 

The bottom panels of Fig. [7J show the corresponding 
average mass accretion histories and compare them with the 
mean predicted MAHs shown in the bottom panels of Fig. [4] 
The latter are the NFW MAHs that result from the (p-2)- 
p C rit(z-2) correlation shown in Fig. [3] The residuals from 
this predicted MAH are clearly similar in shape to those in 
the top panels: in other words, on average, halos whose mass 
profiles deviate from NFW have mass accretion histories that 
deviate from the NFW shape in a similar way. 

Quantitatively, this implies that the best-fit Einasto pa- 
rameters of both MAHs and mass profiles must be corre- 
lated. Since we expect the correlations to be weak (see, e.g., 
Fig. [5} we group halos by mass, concentration, and shape 
parameter (as measured from their mass profiles) before fit- 
ting Einasto profiles to their corresponding average mass 
accretion histories. To prevent possible biases induced by 
numerical-resolution effects the grouping is such that we re- 
tain only well-resolved halos with similar numbers of par- 
ticles, 25, 000 < A200 < 50, 000. Statistical fluctuations are 
reduced by averaging over groups of at least 25 halos, using 
a grid in the c-a plane with a mesh of width 8 log c = 0.079 
and 8a = 0.026. 



The results are shown by the heavy symbols in Fig. [8] 
The left panel shows that the MAH shape parameter is 
clearly correlated with the shape parameter of the mass pro- 
file. Symbols of different colors are used for halos in each of 
the three mass bins, which correspond to different MS. Note 
that they all delineate the same trend, despite the fact that 
they span a range of roughly four decades in virial mass. 
Note as well that parameters corresponding to individual 
halos (shown as grey dots in the figure) correlate less well, 
and that the scatter is larger for the MS-XXL halos. This 
is because individual MAHs are often quite complex, espe- 
cially when major mergers are involved or halos are recently 
assembled and still unrelaxed (even though they pass the re- 
laxation criteria s et out in Sec.l2T2"jl. a s is the case for many 
MS-XXL systems (|Ludlow et al.ll2012l ). These MAHs cannot 
be well approximated by the Einasto shape, thus hindering 
the interpretation of their fit parameters. We therefore focus 
the discussion on the parameters fit to the averaged profiles, 
shown with heavy symbols in Fig. [8] 

The relation between shape parameters is quite weak 
(the 1:1 relation is indicated by the dotted line), implying 
that large variations in MAH shape map onto a narrower 
range of mass profile shapes. As a result, even halos that 
assemble early and over a very short period of time, such as 
those whose MAH is characterized by a ~ 0.5 (see Fig. [6}, 
end up with a ~ 0.25 mass profiles that differ only slightly 
from NFW. This is consistent with earlier findings that halos 
assembled monolithically and without protracted accretion, 
such as those formed in hot dark matter universe s, are never- 
theless well approximated b y NFW profiles (e.g., lHuss et al.l 
ll999l : IWang fc Whitdl2009l) . 
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It also explains why the NFW profile fits rather well ha- 
los formed in hierarchical scenarios other than LCDM. For 
example, the accretion histories of halos formed in scale- free 
scenarios characterized by a power-law spectrum of density 
fluctuations, P(k) oc k n , depends on n, but rather weakly. 
For given n, the MAH shapes are also, on average, inde- 
pendent of halo mas s. Fitting E i nasto profiles to accretion 
histories taken from [Zhao et al.l (|2009f ). we find that halos 
that form from white-noise spectra (n = 0) have MAHs well 
described by a ~ 0.1. For n — —2, on the other hand, the av- 
erage MAH shape is roughly a ~ 0.2. These different MAHs 
result in only a subtle change in mass profile (see left panel 
of Fig. [8]), which would have been undetectable at the nu- 
merical resolution pro bed by earlier work. It m ay, however, 
be behind the claim bv lKnollmann et al.1 HHI) that the in- 
ner slope of the density profile varies systematically with the 
spectral index n. 

The right-hand panel of Fig.[S]shows, on the other hand, 
the correlation between the best-fit Einasto concentration 
parameters of the MAH and mass profiles, for the same set 
of halos shown in the left-hand panel of the same figure. 
This is analogous to Fig. [5] but for Einasto, rather than 
NFW, concentrations. (The NFW c-c correlation is shown 
by a dashed line.) Because of the extra parameter, the re- 
lation between Einasto concentrations depends on a, and is 
indicated by the colored lines in the figure for three values 
of a. The symbol types in both panels are the same, but are 
colored by a in the right-hand panel (see color bar inset). 
Note that, for given a, the c-c relation for MS halos follows 
closely the expected correlations. 

As in the left panel, grey dots in the right-hand panel 
of Fig. [8] correspond to fits to individual halo MAH and 
mass profiles. These clearly follow the same trend as the 
averaged profiles but with larger scatter. In particular, com- 
plexities in the MAH caused by major mergers result at 
times in extreme values for the "concentration" measured 
from accretion histories. These affect in particular large mass 
halos that have rec ently been assembled. As discussed by 
iLudlow et all (|2012l ). many such halos pass the relaxation 
criteria but are actually out of equilibrium and in a par- 
ticular phase of their virialization process. These halos are 
actually responsible for most of the outlier points in Fig. [8] 

We conclude that there is strong evidence for a link 
between the concentration and shape of the mass profile of 
halos and their accretion histories. The correlations are well- 
defined but weak, in the sense that even MAHs whose shapes 
deviate substantially from the mean lead to halos that de- 
part only subtly from the average, NFW-like mass profile. 
This is probably due to the virialization process, where non- 
linear effects lead to a substantial but incomplete erasure of 
memory of the initial conditions from the equilibrium struc- 
ture of a halo. The structural similarity of CDM halos thus 
seems to arise from the mass-independence of MAH shapes 
aided by the homogenizing effect of halo virialization. 



5 SUMMARY AND CONCLUSIONS 

We have examined the mass profile and accretion histories 
of equilibrium cold dark matter halos identified in the Mil- 
lennium Simulation series. As reported in earlier work, halo 
mass profiles are well approximated by NFW profiles which, 



at given virial mass, are characterized by a single parameter, 
the concentration, c = r2oo/f-2- Although in general devia- 
tions from the NFW profile are small, improved fits may be 
achieved using Einasto profiles characterized, at given virial 
mass, by the concentration and an extra shape parameter, 
a. 

Our main finding is that these parameters are strongly 
linked to the accretion history of a halo. The mean density 
within the scale radius, r_2, is directly proportional to the 
critical density of the Universe at the time when the main 
progenitor's mass equals that within r_2. Scaled to these 
characteristic values of mass and density the shape of the 
mass accretion history, expressed as M(p CI i t (z)), is, on av- 
erage, independent of halo mass. Furthermore, this shape 
is nearly identical to that of the enclosed mass-mean inner 
density profile (M((p))) of the halo at z — 0, which can be 
well approximated by the NFW profile. 

This result suggests that the structural similarity of ha- 
los of different mass is related to the fact that their accretion 
histories are independent of mass. It also clarifies how the ac- 
cretion history determines the concentration of a halo; since 
accretion history and mass profile follow the same NFW 
shape, there is a unique correspondence between the "con- 
centration" parameters of either one, as shown in Fig. [5] 

This conclusion is strengthened by the finding that ha- 
los whose mass profiles deviate from NFW and are better 
approximated by Einasto profiles, have mass accretion histo- 
ries that deviate from the NFW shape in a similar way. This 
suggests that the extra shape parameter of the Einasto pro- 
file arises because some halos have accretion history shapes 
that differ from NFW. Indeed, fitting Einasto profiles to 
both M({p)) and M(p cr i t (z)) yields correlated concentra- 
tion and shape parameters. The correlations are clear but 
weak, implying that only halos whose accretion histories de- 
viate strongly from the NFW shape would have mass profiles 
that deviate noticeably from the average, NFW-like shape. 
We ascribe this result to the convergent effects of virializa- 
tion, which partially erase the memory of initial conditions 
from the halo structure. 

Our results suggest that the density profiles of halos 
formed in hierarchical scenarios other than CDM or mono- 
lithically, as in a hot dark matter Universe, are not truly self- 
similar. The deviations from similarity, however, are subtle, 
and a dedicated program of high-resolution numerical sim- 
ulations is needed in order to validate this prediction. 

We may use these findings to predict the dependence of 
halo concentration on mass and redshift, as well as the in- 
fluence of varying the spectral index or the cosmological pa- 
rameters, provided that realistic accretion histories are avail- 
able, either through direct numerical si mulation or through 
well-tested semi-analytic modeling (e.g . , iMonaco et al 1120021 ; 
Ivan den Bosch] |2002| ; Izhao et al.ll2009l ). Such studies would 
help to reveal any shortcomings in our interpretation and 
should shed further light onto the mechanisms responsible 
for CDM halo structure. We plan to address these issues in 
a forthcoming paper. 

Finally, our findings provide some endorsement to the 
many previous studies that have sought a link between 
the final structure of a halo and its evolut i onary history 



(e.g- lBullock et al. ll200ll;IWechsler et alj2003 :l Alvarez et al.l 
[ 20031; IZhaoet al.ll2003l;lTasitsiomi et al.ll2004l ; iLu et al.ll2006l ; 
Izhao et al.ll2009l ; IWong fc Tavlorll2012f ) but still fall short of 
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providing a full account of what determines the structure 
of dark matter halos. Hopefully, the link with accretion his- 
tory we describe here will help to guide future theoretical 
work in order to unravel the mechanism at the root of the 
remarkable structural similarity of dark matter halos. 
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